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CHAPTER  1 .  INTRODUCTION 

This  study  is  to  determine  the  optimal  scheduled  maintenance  policy 
that  performing  a  preventive  replacement  once  the  critical  item  has  reached 
a  specified  age  in  addition  to  failure  replacement  when  necessary.  The 
replacement  action  returns  the  item  to  the  as  new  condition,  thus  con- 
tinuing to  provide  the  same  service  as  the  item  just  been  replaced.  This 
assumption  implies  that  various  costs,  failure  distribution  etc.  used  in 
the  analysis  are  always  the  same. 

A  variety  of  criteria  have  been  used  for  finding  the  optimal  replacement 
age.  In  the  conventionally  approaches,  usually  only  one  criterion  was  con- 
sidered. However,  in  a  real  life,  the  decision  maker  always  wants  to 
attain  more  than  one  objective  or  goal  in  selecting  the  course  of  action. 
Then,  the  problem  becomes  how  to  determine  the  preventive  replacement  age 
of  a  critical  item  subject  to  multiple  criteria. 

The  objectives  of  this  study,  therefore  are:  (i)  to  formulate  the 
model  of  a  replacement  problem  subject  to  multiple  criteria.  (11)  to  solve 
this  problem  by  using  several  multiple  objective  decision  making  methods. 

In  Chapter  2,  the  concepts  of  planned  replacement  are  introduced. 
Among  the  topics  discussed  are:  the  classification  of  maintenance  models, 
the  conditions  required  to  make  preventive  replacement  worthwhile,  factors 
influence  the  preventive  replacement  age,  the  needs  of  multiple  objective 
decision  making. 

In  Chapter  3,  taking  aircraft  as  an  example,  the  problem  is  to 
determine  the  preventive  replacement  age  of  the  critical  item  (e.g.  an 


engine).  Mathematical  models  have  been  developed  for  four  criteria: 
(1)  the  replacement  cost  per  unit  time,  (2)  the  availability  of  the  critical 
item,  (3)  the  mission  reliability,  and  (4)  the  expected  cost  of  failure 
during  the  mission. 

In  Chapter  4,  solutions  are  found  by  using  four  methods  for  multiple 
objective  decision  making:  (1)  the  stricted  selection  method,  (2)  the 
lexicographic  method,  (3)  Waltz's  lexicographic  method,  and  (4)  the  se- 
quential multiple  objective  problem  solving  technique  (SEMOPS).  The  de- 
cision process,  the  implications,  the  advantages  and  disadvantages  of  each 
method  are  discussed  in  the  chapter. 

In  Chapter  5,  the  present  study  is  briefly  summarized.  The  extensions 
of  the  present  study  to  the  multiple  criteria  problem  with  multiple  decision 
variables  and  to  the  group  replacement  problem  are  proposed. 


CHAPTER  2.  CONCEPTS  OF  PLANNED  REPLACEMENT 

2.1  Planned  replacement 

According  to  the  surveys  made  by  Jorgenson  and  McCall  [9]  and  McCall 
[14],  models  of  stochastically  failing  equipment  fall  into  two  classes. 
In  'preparedness'  models  the  equipment  is  subject  to  random  failures  and  the 
actual  state  of  the  equipment  -  good  or  failed  -  is  known  only  with  some 
uncertainty.  The  corresponding  physical  situation  is  that  equipment  is 
required  for  an  emergency,  but  the  equipment  is  not  operated  until  the 
emergency  occurs.  In  'preventive  maintenance'  models,  the  equipment  is 
subject  to  random  failures,  but  the  equipment  is  operated  continuously  so 
that  the  actual  state  of  the  equipment  is  always  known  with  certainty. 

Under  different  knowledge  about  the  distribution  of  times  to  failure 
(the  failure  distribution  is  known,  the  form  of  the  distribution  is  known 
but  the  failure  rate  is  unknown,  and  the  failure  distribution  is  unknown 
etc.)  the  models  can  be  further  classified  by  the  maintenance  procedures 
implemented  into  periodic  policy,  sequential  policy  and  opportunistic  policy. 

The  periodic  policy  in  'preventive  maintenance'  models  with  failure 
distribution  known  has  been  studied  extensively  [1,9].  It  is  noted  that 
the  maintenance  actions  such  as  repair  and  overhaul  can  be  considered  to 
be  equivalent  to  replacement  provided  that  it  is  reasonable  to  assume  that 
such  actions  also  return  the  equipments  to  the  as  new  condition. 

Since  the  earliest  study  on  the  planned-replacement  problem  was  made 
by  Campbell  [5]  in  1941,  many  have  published  on  this  topic.  The  basic 
problems  which  were  discussed  and  analyzed  are:  (1)  what  is  a  valid  and 
good  criterion  to  find  the  optimal  replacement  age?  (2)  when  is  the 
planned  replacement? 


Several  investigators  have  demonstrated  [1,4,6,16]  (Some  with  more 
and  some  with  less  mathematical  rigor)  that  when  a  measure  is  selected  as 
a  criterion,  two  conditions  are  required  to  make  planned  replacement  po- 
tentially worthwhile.  The  first  is  that  the  preventive  replacement  of  a 
item  must  'cost'  less,  in  some  sense,  than  an  failure  replacement.  The 
second  condition  is  that  the  failure  characteristic  of  the  item  must  display 
'wearout';  i.e.  the  failure  rate  must  increase  with  age. 

2.2  Factors  influence  the  replacement  age 

The  following  two  factors  make  the  replacement  models  different; 

(1)  the  measure  selected  as  a  criterion,  and 

(2)  the  failure  characteristic  of  the  item. 

2.2.1  Criterion  for  assessing  the  replacement  age 

The  following  typical  criteria  have  been  used  for  finding  the  optimal 
replacement  age. 

(1)  Expected  cost  per  unit  of  time 

The  cost  of  a  preventive  replacement  and  the  cost 
of  a  failure  replacement  are  considered.  The  objective 
is  to  determine  the  optimal  replacement  age  of  the  equip- 
ment to  minimize  the  total  expected  replacement  cost  per 
unit  time. 

(2)  Maintenance  downtime 

Barlow,  Hunter,  and  Proschan  [3]  first  considered  maintenance 
downtime,  and  derived  an  integral   equation  leading  to  the 
planned  maintenance  interval  which  minimizes  this  time. 


(3)  Availability 

Availability  is  defined  [13]  as  the  fraction  of  the  total 
desired  operating  time  that  the  system  or  item  is  actually 
operable.  It  is  the  principal  measure  of  the  effectiveness 
of  maintained  systems.  The  objective  is  to  determine  the  re- 
placement age  of  the  equipment  to  maximize  the  availability  or 
to  maintain  above  a  specific  level. 

(4)  Mission  reliability 

Mission  reliability  is  the  probability  that  the  system  or  item 
will  operate  properly  without  a  failure  for  a  given  interval 
of  mission  duration,  given  that  it  was  opererable  at  the  begin- 
ning of  the  mission.     In  other  words,  it  represents  the  conditional 
probability  of  nonfailure  of  the  system  for  a  period  of  time  re- 
quired to  complete  a  mission.     Ladany  and  Aharoni   [12]  consider 
it  as  a  criterion  to  determine  the  optimal   replacement  age  of 
a  item  that  there  is  a  requirement  for  a  standard  minimal   value 
of  mission  reliability  to  be  fulfilled  by  the  item  in  a  mission 
of  specific  length. 

(5)  Average  effective  earning 

In  conditions  of  industrial  production  each  equipment  needs  to 
earn  a  revenue.  The  net  revenue  rate  is  measured  as  the  value 
of  output  per  time  unit  less  the  cost  of  producing  that  output, 
exclusive  of  maintenance  cost.  Kay  [10]  maximizes  the  average 
effective  earning  to  determine  the  optimal  schedule  period. 


(6)  Expected  maximum  failure  cost  during  the  mission 

Some  complex  systems  are  so  expensive  or  a  mission  is  so  impor- 
tant or  a  system  failure  will  involve  so  many  lives  that  we 
cannot  afford  to  have  a  major  failure  during  a  mission.     Although 
the  probability  of  occurring  a  major  accident  is  very  low,  min- 
imizing the  expected  failure  cost  during  the  mission  or 
being  below  a  specific  tolerable  value  is  usually  desired. 
Ladany  and  Aharoni   [12]  defined  the  expected  inflight  cost  for 
an  item  with  specific  mission  length  is  the  product  of  the  in- 
flight failure  cost  (including  loss  of  lives,  damage  to  the  whole 
aircraft,  etc.)  and  the  probability  that  the  item  will   fail 
during  the  mission. 

2.2.2     Failure  characteristic 

One  of  the  prerequisites  for  a  planned-replacement  program  is  that 
the  failure  distribution  has  a  wearout  characteristic,  or  equivalently, 
that  the  failure  rate  increases  with  age. 

A  bulk  of  available  failure  analysis  indicates  that  the  exponential 
failure  distribution  is  by  far  the  most  common  of  all   distributions.     Con- 
sequently, we  know  that  a  large  class  of  devices  is  unsuitable  for  planned 
replacement. 

However,  many  researches  [4,11,17]  disclosed  that  a  significant 
fraction  of  electrical   and  mechanical   parts,  especially  that  of  commercial 
or  military  aircraft,  displayed  a  normal,   logarithmic-normal  or  Weibull 
types  of  failure  distribution.     With  increased  discrepancy  of  actual   failure 


distribution  from  the  exponential   failure  distribution  and  using  mainten- 
ance policy  based  on  exponential   failure  distribution,  the  results  of  these 
kind  of  parts  failure  become  more  critical.     Therefore,  in  this  study, 
Weibull   distribution  will  be  used  in  constructing  replacement  models. 

2.3     The  need  of  multiple  objective  decision  making 

In  the  previous  section,  several   typical   criteria  have  been  described. 
Although  they  are  not  exhaustive,  they  demonstrate  that  there  are  many 
criteria  that  can  be  used  to  determine  the  replacement  age.     Traditionally, 
only  one  criterion  has  been  considered.     The  optimal   replacement  age  for 
one  criterion  is  surely  not  an  optimal  solution  for  other  criterion.     At 
most  of  the  cases,   the  solution  is  even  not  a  satisfactory  solution  if  other 
criteria  were  also  considered.     However,   in  a  real   life,   the  decision  maker 
(DM)  always  wants  to  attain  more  than  one  objective  or  goal    in  selecting  the 
course  of  action.     One  of  the  major  reasons  for  the  scarcity  of  multiple 
objective  formulation  and  consideration  in  literature  is  that  until   recently, 
almost  all   the  solution  strategies  developed  involved  a  single  objective 
function. 

The  growing  tendency  to  incorporate  more  and  diverse  criteria  in  a 
system  design  and  the  resulting  difficulty  in  consolidating  such  criteria 
have  provided  an  impetus  to  the  development  of  multiple  criteria  methodology 
in  various  forms. 

In  recent  years,  multiple  objective  analysis  has  been  applied  to  a  wide 
variety  of  problems   [7]  including  water  resources  management,  econometrics 
and  development  planning,  academic  planning,  capital   budgeting,  financial 


planning,   land  use  planning,  manpower  planning,  media  planning,   public 
administration,   transportation  planning,  .and  many  others.     However,  no 
paper  has  been  found  applying  multiple-objective  decision  making  methods 
in  the  field  of  reliability  analysis  and  maintenance  policy  making. 

In  this  study,  a  replacement  age  problem  subject  to  multiple  criteria 
is  developed  and  four  different  kind  of  multiple  objective  decision  making 
methods  are  introduced  and  applied  to  solve  the  problem. 


CHAPTER  3.  STATEMENT  OF  THE  PROBLEM  AND  MODEL  DEVELOPMENT 

3.1  Statement  of  the  problem 

For  every  piece  of  military  or  commercial  equipment,  there  exists 
critical  items  whose  failure  could  result  in  the  shutdown  of  the  whole 
equipment  or  the  hazardous  conditions  for  individuals  using  the  equipment 
or  the  unaccomplishment  of  the  tactical  functions  of  the  equipments.  As 
a  consequence,  it  is  desirable  to  use  multiple  criteria  in  determining  the 
replacement  age  of  such  critical  items.  The  following  four  criteria  are 
most  likely  to  be  considered: 

(1)  the  replacement  cost  per  unit  of  time 

(2)  the  availability  of  the  equipment 

(3)  the  reliability  of  mission  accomplishment 

(4)  the  expected  cost  of  mission  failure 

Reliance  on  a  single  criterion  might  cause  heavy  damage  or  high  opportunity 
cost. 

In  determining  the  maintenance  policy  of  the  critical  item  of  a  air- 
craft, Ladany  and  Aharoni  [12]  selected  the  strictest  result  (i.e.  the 
lowest  replacement  age)  as  the  optimal  replacement  age  by  considering  three 
criteria:  (1)  the  replacement  cost  per  unit  of  time,  (2)  the  mission  reli- 
ability of  the  equipment,  and  (3)  the  inflight  failure  cost.  However,  the 
following  deficiencies  exists  in  their  model  and  the  decision  making  method: 

(1)  The  principal  measure  of  the  effectiveness  of  a  maintained  system, 
availability,  is  not  considered  as  a  criterion  in  the  model. 

(2)  The  equation  of  the  replacement  cost  per  unit  of  time  did  not 
take  account  of  replacement  times  of  the  preventive  and  the 
failure  replacements. 
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(3)     Using  the  strictest  selection  method,  the  decision  maker  is  es- 
sentially detached  from  the  decision  process. 

In  this  study,  in  addition  to  the  strictest  selection  method,  several 
methods  for  multiple  objective  decision  making  are  used  for  obtaining  the 
optimal   scheduled  maintenance  policy  based  on  four  criteria. 

3.2     Model   development 

The  problem  is  to  determine  the  preventive  replacement  age  for  a 
critical   item  of  a  complex  system  (e.g.  an  aircraft)  based  on  the  following 
multiple  criteria: 

(1)  the  replacement  cost  per  unit  time 

(2)  the  mission  reliability  of  a  complex  system 

(3)  the  expected  mission  failure  cost 

(4)  the  availability  of  the  item 

Notations: 


A(t  )  =  the  availability  of  the  critical  item 

C-,    =  the  cost  of  a  preventive  replacement 

Co    *  the  cost  of  a  failure  replacement 

C,    =  the  cost  of  a  mission  failure  (including  loss  of  lives,  damage 

to  whole  system  etc.) 
C(t  )  =  the  total  expected  replacement  cost  per  unit  of  time 
Ef    =  the  expected  mission  failure  cost 
f(X)   =  the  failure  density  function  of  the  critical  item 
F(X)   =  the  unreliability,  F(X)  =  1  -  R(X) 
F(X,H)  =  the  probability  of  failure  between  X  and  X+H,  given  it  is 

not  failed  at  X 


n 


H  =  the  mission  length 

CO 

R(X)       =  the  reliability,  R(X)  =  /  f (X)dt  . 

X 

R(X,H)  =  the  mission  reliability,  R(X,H)  =  1  -  F(X,H) 

t-,  =  the  mean  replacement  time  for  the  preventive  replacement 

t2  =  the  mean  replacement  time  for  the  failure  replacement 

t  =  the  preventive  replacement  age 

Assuming  that  the  planning  horizon  is  very  long,  and  that  an  item 
which  is  replaced  behaves  like  a  new  item  with  zero  length  of  life,  the 
above  mentioned  four  criteria  are  evaluated  as  follow: 
(1)  Replacement  cost  per  unit  of  time 

The  replacement  policy  is  to  perform  a  preventive  replacement  once  the 
critical  item  has  reached  a  specified  age  t  ,  and  a  failure  replacement 
when  it  is  necessary.  This  policy  is  illustrated  in  Fig.  3.1. 

There  are  two  possible  cycles  of  operation: 
One  cycle  being  determined  by  the  critical  item  reaching  its  preventive 
replacement  age  t  ,  the  other  being  determined  by  the  failure  occuring 
before  the  preventive  replacement  age.  These  two  possible  cycles  are 
illustrated  in  Fig.  3.2.  The  total  expected  replacement  cost  per  unit  time, 
C(tp),  is: 

„,  x  .  total  expected  replacement  cost  per  cycle 
^  p'  expected  cycle  length 

Total  expected  replacement  cost  per  cycle 
=  cost  of  a  preventive  replacement  cycle  •  probability  of  a  preventive 
replacement  cycle  +  cost  of  a  failure  replacement  cycle  •  probability 
of  a  failure  replacement  cycle 
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Fig.  3.2  Preventive  and  failure  replacement  cycle 
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=  C1  R(tp)  +  C2  F(tp) 

Expected  cycle  length 

=  length  of  a  preventive  replacement  cycle  .  probability  of  a  pre- 
ventive replacement  cycle  +  expected  length  of  a  failure  replacement 
cycle  .  probability  of  a  failure  replacement  cycle 

=  (t  +  t,)R(t  )  +  (expected  length  of  a  failure  replacement  cycle)  F(tp) 

If  a  preventive  replacement  occurs  at  time  t  ,  then  the  mean  time  to 
failure  is  the  mean  of  the  shaded  portion  of  Fig.  3.3.  Since  the  unshaded 
area  is  an  impossible  region  for  failure,  the  mean  of  the  shaded  area  is: 

/  p  tf(t)dt 


1  -  R(tp) 


Then: 


the  expected  length  of  a  failure  replacement  cycle 
t  tn 

/  p  tf(t)dt  /  p  tf(t)dt 

~   i  -  R(tj    +  h  -  ^TG  h 


the  expected  cycle  length 


=  (tn  +  t,)  R(tn)  + 


/  p  tf(t)dt 
0  +  t 

~KtJ  +tJ 


F(tp) 


tpR(tp)  +  trR(tp)  +  /  p  tf(t)dt  +  t2F(tp) 


0 


Since 


t  t 

/    R(x)dx  =  t-R(t)  +  /     tf(t)dt 
0  0 
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Service 
Time 


Fig.  3.3.  Mean  time  to  failure 
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t 

The  expected  cycle  length  =  /  p  R(t)dt  +  t1R(tp)  +  t2F(tp) 

Then 

t     W+c^y (3]) 

/  p  R(t)dt  +  t^Ctp)  +  t2F(tp) 

This  relates  the  preventive  replacement  age,  t  ,  to  the  total  expected  re- 
placement cost  per  unit  time. 

(2)  Availability  of  the  critical  item 

Availability  is  a  principal  measure  of  the  effectiveness  of  a  main- 
tained system.  Availability  is  defined  [13]  as  the  fraction  of  the  total 
desired  operating  time  that  the  system  is  actually  operable,  or  it  can  be 
defined  as  the  ratio  of  uptime  to  total  time. 

When  the  service  time  of  the  critical  item  is  considered,  the  avail- 
ability of  the  item,  A(t  ),  is: 

./.  ,  _  expected  operating  length  in  a  cycle 
^  p'      expected  cycle  length 

The  expected  operating  length  in  a  cycle 

=  operating  length  in  a  preventive  cycle  •  probability  of  a  pre- 
ventive cycle  +  operating  length  in  a  failure  cycle  •  probability 

of  a  failure  cycle 

t 

=  t_R(t)  +  /  p  tf(t)dt 
P   P    0 

t 
=  J  p  R(t)dt 
0 
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as  the  expected  cycle  length 

t 
r 

o 
The  availability  is 


•  /  P  R(t)dt  +  t1R(tp)  +  t2F(tp) 


tn 

/  p  R(t)dt 

A(tp)   =  -j °- (3.2) 

/  p  R(t)dt  +  t,R(t_)  +  t2F(t) 

This  is  a  measure  of  the  critical   item's  availability  because  it  gives 
the  probability  the  system  will  be  available  for  operation  if  it  is  in  an 
ideal   support  environment  without  consideration  for  downtime  other  than 
t|   and  t2»  such  as  administrative  time,  supply  time  etc. 

(3)     Mission  reliability 

A  mission  is  an  interval   of  activity  of  length  H  which  starts  at  age  X 
of  the  item  and  ends  at  age  X+H.     During  the  interval   of  length  H,   the 
equipment,  say  aircraft,  may  even  perform  more  than  one  flight,  but  all 
flights  belong  to  the  same  mission  (e.g.,  a  commercial   flight  with  landings 
for  refuelling) . 

The  conditional   probability  of  a  failure  during  a  mission  of  length  H 
of  a  critical   item  given  it  is  operable  at  age  X,   F(X,H),  is  the  probability 
that  it  will   fail   between  X  and  X+H,  divided  by  the  probability  it  does  not 

fail   until   X,   i.e., 

X+H 

/        f(t)dt 

c/Y  u\       0 _  F(X+H)   -   F(X)  _  1-R(X+H)    -  1+R(X) 

F(X,H) ^  =  Rfx)  =  Rfe 

_  ,        R(X+H) 

~  ]  ""RTxT 


18 


Then,  the  mission  reliability,  which  is  the  probability  for  the  item  to 
succeed  during  the  mission  of  length  H,  if  it  started  at  age  X  is 

R(X,H)  =  l-F(X.H)  =^yi  (3-3) 

Reliability  is  a  monotonic  decreasing  function  of  time.  The 
mission  reliability,  R(X,H),  will  be  the  minimum  at  the  preventive  replace- 
ment age,  that  is  R(t  ,H). 

(4)  Expected  mission  failure  cost 

When  a  system  is  very  expensive,  or  a  mission  is  very  important,  or  a 
system  failure  will  involve  many  lives,  it  cannot  afford  to  have  a  major 
failure  during  a  mission.  The  cost  of  a  mission  failure  includes  the  loss 
of  lives,  the  damage  to  whole  system  etc.  The  expected  mission  failure  cost 
during  a  mission  of  length  H,  Ef,  is  the  product  of  the  mission  failure  cost 
and  the  probability  that  the  item  starting  the  mission  at  age  X  will  fail 
during  the  mission;  i.e. 

Ef  -  C3  •  F(X,H)  •  Cj[l  -  *$})  (3-4) 

Unreliability  is  a  monotonic  increasing  function  of  time.  The  mission 
unreliability,  F(X,H),  will  be  the  maximum  at  the  preventive  age,  that  is 
F(tp,H). 

3.3  Case  of  Weibull  distribution 

When  the  failure  of  an  item  is  distributed  according  to  the  Weibull 
distribution,  the  failure  density  function,  the  failure  function,  the 
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reliability  function  and  the  hazard  function,  respectively,  as  found  in 
Von  Alven  [18]  are: 


f(t)  =  f  t8"1  exp(-ItB)  (3.5) 

F(t)  =  1  -  exp(-ltB)  (3.6) 

R(t)  =  exp(-ltB)  (3.7) 


and 


Z(t)=|tB-'  (3.8) 


where  A  and  B  are  parameters  of  the  distribution  which  has  a  mean  of 

1 
E(t)  =  AB  r  (1+  1) 

and  a  variance  of 

2 


V(t)  =  AB  {r(l  +1)  =  [r(l+  1)}; 


Substituting  eqs.  (3.6),  (3.7)  into  eq.  (3.1),  we  obtain  the  replacement 
cost  per  unit  of  time  as 

C(tp)  .  -, ^xP(-}tBp)+C2[l-exp(-ltB)]  ^  ^ 

/  p  exp(-  I  tB)dtp  +  t]   exp(-  J-  tB)  +  t2[l-exp(l    -  1  tB)] 

Substituting  eqs.   (3.6)   and  (3.7)   into  eq.    (3.2)   for  the  availability  model, 
we  get 
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A  exp(-  1  tB)dtp 


A(tJ  = * (3.10) 

/ 
0 


/Up  exp(-  \  tBp)dtp  +  t,  exp(-  \  tB)  +  t2[l  -  exp(-  }  tB)] 


Substituting  eq.   (3-7)  into  eq.   (3.3)  for  the  case  of  the  mission  reliability 
model ,  we  get 

exp[-  I  (t_  +  H)B] 
R(tp,H)  =  ^— J2 (3.11) 

exp  (-  I  tBp) 

Substituting  eq.  (3.7)  into  eq.  (3.4)  for  the  case  of  the  expected  mission 
diss  aster  cost  case,  we  get 

r        exp  [-  I  (t  ♦  H)B]  s 
E   =  C3  •  F(tp  H)  =  C3{l LJ- }  (3.12) 

Numerical  example 

A  realistic  numerical  example  is  taken  from  Bell,  Kamins  &  McCall  [3]. 
An  aircraft  engine  is  found  to  have  failure  characteristics  which  are 

Q 

closely  represented  by  a  Weilbull   distribution  with  A  =  2.6954  x  10  , 
B  =  3.0. 

It  is  also  reasonable  to  assume  the  following  numerical   values: 

C1   =  $25,000 

C2  =  $37,500 

C3  =  $2,500,000 

t1   =  8  hrs 

t2  =  16  hrs 

H     =  16  hrs 
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The  figures  of  failure  density  function,  the  failure  function,  the  relv 
ability  function  and  the  hazard  function  are  shown  in  Figs.  3.4,  3.5, 
3.6  and    3.7,  respectively. 
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CHAPTER  4.  REPLACEMENT  AGE  OF  THE  CRITICAL  ITEM 
SUBJECT  TO  MULTIPLE  CRITERIA 

In  this  chapter,  the  optimum  replacement  ages  of  a  critical  item  for 
an  aircraft  subject  to  multiple  criteria  are  obtained  by  several  methods. 
One  method  (the  strictest  selection  method)  was  proposed  by  Ladany  and 
Aharoni  [12].  The  other  newly  developed  methods  for  multiple  objective 
decision  making  are  applied  to  this  field  for  the  first  time. 

The  numerical  example  shown  in  section  3.3  is  repeatly  used  to  demon- 
strate the  decision  making  methods.  The  results  obtained  from  these 
methods  are  compared  and  discussed. 

4.1  The  strictest  selection  method 

With  multiple  criteria,  it  will  provide  the  critical  item  with  dif- 
ferent preventive  replacement  ages  for  each  criterion.  Ladany  and  Aharoni 
[12]  proposed  that  the  optimal ity  of  maintaining  an  aircraft  requires  the 
application  of  the  strictest  maintenance  policy,  i.e.,  the  criterion  of 
the  shortest  preventive  replacement  age  will  be  applied. 

To  apply  the  strictest  selection  method,  the  optimal  preventive  re- 
placement age  for  each  criterion  has  to  be  found  out  first.  By  the  models 
developed  in  chapter  3,  the  optimal  preventive  replacement  ages  for  the 
four  criteria  are  found  as  follow: 
(1)  Replacement  cost  per  unit  of  time 

The  minimal  replacement  cost  per  unit  of  time  is  obtained  by  equalling 
the  first  derivative  of  eq  (3.1)  to  zero.  The  resulting  condition  for 
the  optimal  preventive  replacement  time,  t*,  is 
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C„t,-C,t, 


Z(t*)   /  P  R(t)dt  -   F(t*)  =  p-4-  -  Z(t*)   [  V-C1   ^ 


(4.1) 


where  Z(t  )  is  the  hazard  function.     In  the  case  of  Weibull   distribution, 
t*  is  obtained  when  the  following  equation  is  fulfilled. 


i  t*  B-1 


1   ,B, 


/'  exp<-  ±  t°)dtn  + 


C2tl'Clt2 


A  V     P         C2-C1 


-  1   +  exp(-itJB) 


C2"C1 


=  0 


(4.2) 


With  the  numerical   example  given  in  section  3.3,  the  curve  of  replacement 
cost  per  unit  of  time  C(t  )  is  shown  in  Fig.  4.1.     The  optimal   solution 
is  found  by  using  a  numerical   integration  and  the  golden  section  method. 
The  optimal   preventive  replacement  age  of  this  specific  numerical   example 
is  1455  hours. 
(2)     Availability  of  the  critical   item 

Using  availability  as  a  criterion  for  deciding  the  preventive  replace- 
ment age  of  the  item,  either  requires  to  find  a  t     value  that  maximize  the 

r 

availability  or  to  find  a  heighest  t  value  satisfying  a  standard  minimum 
value  A(X)mini  i.e., 


A(X)nHn   C  A(tJ 

'mm  —  x  p' 


In  the  case  of  maximizing  the  availability,  the  optimal  preventive  replace- 
ment time  t*  is  obtained  again  by  equating  the  first  derivative  of  eq.  (3.2) 
to  zero.  The  resulting  condition  for  t*  is 


'p 


Z(tJ)  /  p  R(t)dt  =  j-f--   R(t*) 
v     0  z2     1     p 


(4.3) 
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where  Z(t  )   is  the  hazard  function. 

ForWeibull  distribution,  the  optimal  preventive  replacement  age  t*  is 

obtained  when  the  following  equation  is  fufilled: 

|t^^exp(4tB)dtp-^-exp(-}^,  (4.4) 

With  the  numerical  example  given  in  section  3.3,  the  curve  of  availability 
of  the  item  is  shown  in  Fig.  4.2.  The  optimal  preventive  replacement  age 
of  this  specific  numerical  example  is  1129  hours. 
(3)  Mission  reliability 

It  is  desired  that  the  probability  of  failure  of  the  critical  item 
will  not  exceed  the  maximum  acceptable  risk.  So  as  a  criterion  for  selecting 
the  preventive  replacement  age  of  the  item,  it  usually  requires  the  item  to 
fulfill  a  standard  minimal  value  R(X,H)min  in  a  mission  of  length  H.  Then 
R(X,H)  •  becomes  the  limiting  condition  for  determining  t*  for  the  critical 
item.  It  is  required  that 

R(tS+H)  (A   <0 

R(X,H)  •  <  B (4-5) 

'■In-  R(t*} 

and  the  highest  value  of  t*  fulfilling  about  equation  should  be  chosen. 
For  the  numerical  example  with  Weibull  distribution  given  in  section  3.3, 
the  curve  of  mission  reliability  of  the  item  is  shown  in  Fig.  4.3.  Setting 

R(X,H)    at  0.985,  the  optimal  preventive  replacement  age  of  this  specific 

v   'mm 

numerical  example  is  913  hours. 
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(4)  Expected  inflight  failure  cost 

It  is  obvious  that  if  the  critical  item  is  replaced  with  a  preventive 
replacement  just  before  the  mission,  the  new  item  will  have  an  age  of  zero. 
The  probability  of  failure  during  the  mission  will  be  minimal  so  is  the 
expected  mission  failure  cost.  So  the  condition  for  finding  t*  for  the 
critical  item  is  obtained  from  the  requirement  that  the  expected  mission 
failure  cost  should  equal  the  preventive  replacement  cost  of  the  item  just 
prior  to  the  mission  plus  the  expected  mission  failure  cost  of  newly  re- 
placed critical  item,  i.e., 

C3F(t*,H)  =  C3F(0,H)  +  C1 
that  is 

E =  R(H)  -  p1  C4.6) 

R(t*)        L3 

For  the  numerical  example  with  Wei  bull   distribution  given  in  section  3.3, 
the  curve  of  expected  mission    failure    cost  C3F(t  ,H)   is  shown  in  Fig.  4.4. 
By  eq.   (4.6)  we  can  find  the  optimal   preventive  replacement  age  at  743  hours 
with  minimal  expected  mission    failure  cost  $25,005  which  is  equal   to  the 
cost  that  can  occur  by  replacing  the  critical   item  just  before  the  mission. 
The  following  results  in  hours  for  the  optimal   preventive  replacement 
age,  t*,  were  found  for  the  critical   item  according  to  the  four  criteria: 
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Criteria  *£  (nr-) 


E. 


replacement  cost  per  unit  of  time  1,455 

availability  1,129 

mission  reliability  913 

mission    failure   cost  743 

The  lowest  replacement  age,  743  hours,  has  been  chosen  in  order  to  assure 
that  (a)  at  least  the  minimal  mission  reliability  is  fulfilled,  (b)  the 
expected  mission    failure   cost  is  not  higher  than  the  total   cost  consisting 
of  the  expected  mission    failure   cost  of  a  new  item,  and  the  preventive 
replacement  cost,   (c)  the  lowest  replacement  cost  per  unit  of  time  and  the 
highest  availability  are  obtained,  given  that  (a)  and  (b)  are  fulfilled. 

The  optimal   replacement  age  obtained  by  the  strictest  selection  method 
also  demonstrated  in  Fig.  4.5. 

4.2     Lexicographic  method 

This  method  requires  that  the  objectives  are  ranked  in  order  of  impor- 
tance by  the  decision  maker  (DM).     The  satisfying  solution  obtained  by  this 
method  is  one  which  maximizes   the  objectives  starting  with  the  most  important 
and  going  down  according  to  the  order  of  the  importance  of  the  objectives  [7]. 

Assume  that  a  multiple  objective  decision  problem  consists  of  n  decision 
variables,  m  constraints  and  K  objectives.     Let  the  subscripts  of  the  objec- 
tives indicate  not  only  the  components  of  the  objective  vector,   f_(X_) ,  but 
also  the  priorities  of  the  objectives,   i.e.,   f,(X)   is  the  first  component 
of  f(X)  and  the  most  important  objective,   f2(X)   is  the  second  component  and 
the  second  most  important  objective,  and  so  on.     Then  the  first  problem  to 
be  solved  is: 
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Max 


MX) 


(4.7) 


s.t.    g.(X)  <  0    j  »  1,2,  ...,  m  J 

J 


Let  ft   be  the  solution  to  the  above  problem.  If  this  problem  gives  a  unique 
X^  for  ft,   the  solution  is  considered  as  the  satisfying  solution  to  the  entire 
problem.  Otherwise,  the  second  problem  is  to  be  solved,  i.e., 


Max    f2(X) 


s.t.    g.(X)  <.  0   j  =  1 ,  2,  ...,  m  • 

J 


f ! (X)  "  f { 


(4.8) 


Let  f*  be  the  solution  to  (4.8).  If  (4.8)  gives  a  unique  X  for  f|,  the 
solution  is  the  satisfying  solution  'to  the  entire  problem.  Otherwise  the 
procedure  is  repeated  until  all  K  objectives  have  been  considered.  In 
general,  the  ith  problem  is: 


Max    f  i  ( X) 

s.t.    g. (X_)  <_  0   j  =  1 ,2,  . ..,  m 

J 


f£(X)  =  f J  A  =  1,2,  ...,  1-1 j 


(4.9) 


Since  the  procedure  is  terminated  when  a  unique  solution  is  obtained  at  the 
ith  problem;  the  solution  will  be  the  satisfying  solution  to  the  entire 
problem;  the  objectives  ranked  less  important  than  f^(X.)  are  ignored  by 
this  method.  The  rationale  for  this  method  is  that  individuals  tend  to 
make  decisions  in  this  manner. 
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The  Replacement  Age  Example 

Assume  the  DM  ranks  the  importance  of  four  criteria  as:     (1)  mission 
reliability  (2)  expected  cost  of  inflight  failure  cost  (3)  replacement  cost 
per  unit  of  time  (4)  the  availability  of  the  equipment.     Then  the  first 
problem  to  be  solved  is: 


Find 


s.t. 


R(tp,H)   >R(tpfH)m1n 


S-0 


(4.10) 


where  R(tn,H)    •     =  0.985. 
v  p     'mm 

The  solution  is  t     <_  913  hrs   (see  Fig.  4.6.)     Since  this  is  not  a 

unique  solution,  the  second  problem  to  be  solved  is: 


Find  t 


s.t.  C3F(t  ,H)  <_  C3F(0,H)  +  C] 


R(t  ,H)       ^  0.985 


%     io 


(4.11) 


The  solution  to  (4.11)   is  t     <_  743  (see  Fig.   4.6).     Since  this  is  still 
not  a  unique  solution,  the  third  problem  to  be  solved  is: 
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Max 

{V 

s.t. 


(4.12) 


C(t) 


R(t  ,  H)  >  0.987 


C3F(t   ,H)   <  C3F(0,H)  +  C1 


t     >  0 
P  - 


The  solution  to  (4.12)  is  t    =  743  hrs.     For  this  is  a  unique  solution, 
the  procedure  is  terminated.     The  t    =  743  hr  is  the  satisfying  solution 
to  the  entire  problem.     The  less  important  criterion,  availability  of  the 
equipment,  is    ignored   by  this  method.     The  optimal   replacement  age  of  this 
specific  example  obtained  by  the  lexicographic  method  is  illustrated  in 
Fig.  4.6.     Since  there  are  four  criteria  in  the  problem,  the  DM  may  rank 
their  importance  in  41  =  24  different  ways.     The  results  for  these  24 
different  ways  of  rankings  are  shown  at  Table  4.1.     Results  in  Table  4.1 
indicate  that  for  this  specific  numerical  example,  the  satisfying  solution 
is  dominated  by  the  first  priority  or  the  first  and  second  priorities. 
It  is  noted  that  the  satisfying  solution  will  be  different  if  the  priorities 
of  the  four  criteria  are  changed.     Since  the  solution  is  very  sensitive  to 
the  ranking  of  the  objectives  given  by  the  DM;  the  analyst  should  exercise 
caution  in  applying  this  method  when  some  objectives  are  of  nearly  equal 
importance. 


4.3    Waltz's  lexicographic  method 

A  variation  of  the  method  proposed  by  Waltz  [19]  may  reduce  the  sensi- 
tivity of  the  DM's  priority  of  criteria.     After  the  first  objective  is 
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Table  4.1  The  results  of  different  rankings  of  criterion  importance 
by  lexicographic  method 


Priority  1 

Priority  2 

Priority  3 

Priority  4 

t*(hr) 

c(tp) 

Those  criteria  in 
this  area  will 
not  affect  the 
solution 

1455 

A(tp) 

1129 

R(tp,H) 

A(tp) 

913 

R(tp,H) 

c(tp) 

913 

R(tD,H) 

C3,F(tp5H) 

743 

C3F(tp,H) 

743 
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maximized,  the  second  objective  is  maximized  subjected  to  keeping  the 
first  objective  within  a  certain  percentage  of  its  optimum.  The  third 
objective  is  then  maximized  keeping  the  first  two  within  a  certain  per- 
centage of  the  optimum  values  found  in  the  previous  step.  The  ith  problem 
then  is: 


Max 


Ud) 


s.t.    q.{X)   <  0    j  =  1,2,  ...,  m 


f£(D  -  f*-  \*     I'   1,2,  ...,  j-1 


(4.13) 


where  <5£'s  are  tolerances  determined  by  the  DM. 

The  Replacement  Age  Example: 

Now,  if  the  DM  decides  that  the  importance  of  four  criteria  are  ranked 
as  (1)  the  replacement  cost  per  unit  of  time,   (2)  the  availability  of  the 
equipment,   (3)  the  mission  reliability,  and  (4)  the  expected  inflight 
failure  cost,   then  the  first  solution  to  be  solved  is: 


Min 

{V 


s.t, 


C(tp)  ' 


tp-° 


(4.14) 


The  solution  is  t*  =  1455  hrs,  and  C(t  )  =  28.92  $/hr  (See  Fig.   4.7). 
However,  the  DM  decides  the  replacement  cost  per  unit  time,   C(t  ),   is 
satisfactory  when  it  is  less  than  30.5  $/hr.     Then  the  second  problem  to  be 
solved  is: 
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Max. 

{V 


s.t. 


A<y 


C(t  )  <  30.51-  (4.15) 


S-0 


The  solution  to  (4.15)   is  t*  =  1129  hrs.   and  A(t*)  =  0.9888  (see  Fig.  4.7) 

r  r 

The  DM  thinks  that  it  is  acceptable  when  the  availability  of  the  item  is 
higher  than  0.9875.  Then  the  third  problem  to  be  solved  is: 


ax. 

V 

R(tp.H) 

.t. 

C(t  )  £  30.5 

A(t  )  >  0.9875 

tp-° 

(4.16) 


The  solution  to  (4.16)   is  t*  =  1057  hrs,   R(t*,H)  =  0.98  (See  Fig.   4.7). 
At  this  point,   the  DM  find  that  the  value  of  mission  reliability  can  not 
be  relaxed.     So  a  unique  solution  is  obtained  and  the  procedure  is  termin- 
ated.    The  satisfying  solution  is: 


t*  =  1057  hrs. 


C(t  )  =  30.5  $/hr 
A(t  )   =  0.98877 


R(t  ,H)  =  0.98 


C3F(t  ,H)   =  $50000 
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It  is  noted  that  the  sensitivity  of  the  priority  ranking  to  the  solution 
is  reduced.     The  optimal   replacement  age  of  this  specific  example  obtained 
by  the  Waltz's  lexicographic  method  is  illustrated  in  Fig.  4.7. 

4.4     SEMOPS  -  An  Interactive  Method 

One  class  of  multiple  objective  decision  method  is  generally  referred 
as   'Interactive  Methods'.     It  is  rely  on  the  progressive  definition  of  the 
DM's  preferences  along  with  the  exploration  of  the  criterion  space.     Much 
work  has  been  done  recently  on  this  class  of  methods   [7].     The  progressive 
definition  takes  place  through  a  DM-Analyst  or  DM-Machine  dialogue  at  each 
iteration.     At  each  such  dialogue,  the  DM  is  asked  about  some  trade-off  or 
preference  information  based  upon  the  current  solution  (or  the  set  of 
current  solutions)   in  order  to  determine  a  new  solution.     These  methods 
assume  that  the  DM  is  unable  to  indicate   'a  priori'   preference  information 
due  to  the  complexity  of  the  problems.     However,  the  DM  is  able  to  give  a 
preference  information  on  a  local   level   to  a  particular  solution.     As  the 
solution  process  progresses,  the  DM  not  only  indicates  his  preferences 
but  also  learns  about  the  problem. 

A  Sequential  Mul tiobjective  Problem  Solving  technique  (SEMPOS)  was 
proposed  by  Monarchi,  Kisiel   and  Duckstein  [15].     The  DM  is  allowed  to  trade 
off  one  objective  vs  another  in  an  interactive  manner.     Some  implicit 
trade-off  information  in  the  form  that  the  DM  be  able  to  indicate  whether 
he  is  satisfied  at  the  current  achievement  level   are  generated  by  the 
mechanism  of  the  method.     SEMOPS  cyclically  uses  a  surrogate  objective 
function  based  on  goals  and  the  DM's  aspirations  about  achieving  the  objective. 
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The  goal   levels  are  conditions   imposed  on  the  DM  by  external   forces  and 
the  aspiration  levels  are  the  attainment  levels  of  the  objectives  that  the 
DM  personally  desires  to  achieve.     (Goals  do  not  change,  but  aspiration 
levels  change  as  each  iterative  cycle  goes.) 

Let  AL  =  (AL-, ,   ...,  ALj  be  the  DM's  aspirations  levels,  and 
f(X)  =  {f-,(Xj,   ...,  fJl)}  be  the  multiple  objectige  functions.     First, 
a   'relevant  range'  of  f_(Xj   for  each  goal   is  chosen  as  [f«(X)ii  f-j(X.)i|]- 
The   'relevant  range'   is  the  interval   that  the  DM  thinks  the  f. (X)  will   be 
in  of  the  problem  he  is  considering.     By  transforming  the  original   response 
surface  f.(X)  to  Y-(X)  where  Y^X)  is  defined  on  the  interval   (0,1].     Then 
for  any  goal : 

f-j(X)   -  fj(A)L 
Yi(-}  =  ft(X)u  "  fi(X)L 

Y.j(X_)  is  in  the  interval  (0,1].  In  examining  the  results  of  each  iter- 
ation, it  will  be  necessary  to  confirm  that  f.  (X)  has  indeed  remained  in 
the  specified  interval. 

The  AL_  may  be  transformed  into  A  by  the  same  transformation 

_  *,  -  yx)L 

Ai  "  f1{x)u  -  fj(x)L 

Then  A  is  in  [0,1] 

The  dimensionless  indicator  of  achievement,  d,  is  defined  for  five 
types  of  objectives: 
(1 )  at  most 

f^X)  <  AL1;    di  =  fi(X)/ALi 


(2)  at  least 


f^X)  >  AL^  d.   =  Al^/f^X) 


(3)  equals 


.ifJ^'i® 


'ica-*,i      "i  -  ?  (f^ir +  ?q 


(4)  within  an  interval 


AL 


iLlfiODiAL^;  di   = 


(5)     outside  an  interval 


ALiU         VALiL     ,   V^l 


ALiL  +  ^luJlW      ALiuJ 


46 


or 


*l(X)  lAL.L 


f^X)   >ALiL) 


di   = 


ALiL+ALiU 


AL 


iU 


\    r 

1 

ALa 

fi(A) 

'  ALiU    J 

Jlf,(x) 

Types   (1),   (2)  and  (4)  are  the  most  common.     In  each  instance,    values  of 
di  <_  1   imply  that  the  goal   is  satisfied. 

The  algorithm  generates  information  under  the  guidance  of  the 
DM  so  that  he  can  make  a  decision.     Information  concerning  the  interrela- 
tionships between  objectives  is  in  terms  of  how  achievement  or  non- 
achievement  of  one  objective  affects  the  aspiration  levels  of  other  ob- 
jectives. 

The  mechanism  whereby  information  is  generated  for  the  DM  to  evalu- 
ate is  the  cyclical   optimization  of  a  surrogate  objective  function  s. 
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s  =  ^   d 
teT'  Z 


Where  T'  is  the  subset  of  the  set  of  T  objectives  as  those  objectives 
making  up  s  at  a  given  iteration  of  the  decision -making  process.  It  is 
important  to  note  that  this  optimization  only  provides  information  to  the 
DM;  it  does  not  solve  the  decision  problem.  The  word  surrogate  is  used 
in  recognition  of  the  fact  that  the  true  preference  function  of  the  DM 
is  unknown.  The  value  of  each  d,  in  s  reflects  whether  the  t-th  ob- 
jective has  been  satisfied;  unsatisfied  objectives  have  values  d^  >  1 . 
Of  course,  the  nonlinearity  prevents  direct  comparison  among  the  values  of 
d. 

Operationally,  SEHOPS  is  a  three-step  algorithm  involving  setup, 
iteration  and  termination.  Setup  involves  transforming  the  original 
problem  into  a  principal  surrogate  objective  function  problem  and  a  set 
of  auxiliary  problems  involving  surrogate  objective  functions.  The  iter- 
ation step  is  the  interactive  segment  of  the  algorithm  and  involves  a 
cycling  between  a  optimization  phase  (by  analyst)  and  an  evaluation  phase 
(by  the  DM)  until  a  satisfaction  is  reached,  which  terminates  the  algorithm 

The  first  iteration,  i  =  1,  solves  the  principle  problem,  and  a  set 
of  T  auxiliary  problems  formed  as  follows,  where  the  aspiration  level  of 
each  objective  is  given  as  the  goal  of  each  objective,  i.e.,  AL.  =  b. , 
i  =  1 ,2,  ...,  T. 
The  principle  problem: 

T 

I 
t=l 


min   S-,  =  I    dt' 


s.t.   x  e  X 


(4.17) 
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The  set  of  auxiliary  problem,  i   =  1,2,  ...,  T, 


mm 


'u"  I    dt  1 

\i      t=]     z 


S.t.    x  e  X 


f,(x)  iAL. 


(4.18) 


Solving  (4.17)  and  (4.18)  forms  the  optimization  phase.     The  resulting 
policy  vector  and  objectives  for  the  principal   problem  and  the  set  of 
auxiliary  problems  are  presented  to,  and  are  used  in  the  evaluation 
phase  by  the  DM.     This  information  provides  the  effects  of  changing  the 
aspiration  level   for  any  of  the  goals.     Some  simple  graphical   represen- 
tation of  the  information  can  provide  the  DM  with  rough  estimates  of  the 
interaction  among  the  goals  that  takes  place  through  the  constraint 
set  as  his  aspiration  levels  changes. 

In  general,  the  ith  iteration  solves  the  following  principle 
problem,  and  a  set  of  auxiliary  problems. 
The  principle  problem: 


mm 


s.t. 


'i "  X  dt 

tsl 


fj(x)    >  ALj,        j    e    (J-V)j 


(4.19) 


The  set  of  auxiliary  problems,  i   e  T1  (the  number  of  T'  =  T  -  i  +  1) 
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mm 


s.t. 


U       teT' 


(4.20) 


x  E   X 


fj(X)  lALj  for    Yj,     j  e  (T-T') 


fA(X)  1  AL  ,         for  one  £,         a  e  T' 


The  optimization  phase  solves  (4.19)  and  (4.20).     The  resulting  solutions 
are  used  in  the  evaluation  phase  and  a  guidance  is  given  by  the  DM  for 
the  next  iteration  cycle.     The  search  continues  until  a  satisfactory 
solution  is  found. 

The  Replacement  Age  Example: 

In  this  study  of  finding  the  replacement  age  involves  four  goals 
and  one  non-negative  decision  variable.     Their  various  formulations  are 
expressed  below: 

Goals  and  criterion  functions 


h   =  C(tn)   = 


C-,R(tp)  +  C2F(tp) 


/  P  R(t)dt  +  t^(t  )  +  t2F(t  ) 


<  AL] 


/  p  R(t)dt 


h  -  A<V  =  - 


/  p  R(t)dt  +  t^(t  )  +  t2F(t  ) 


>  AL2 
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R(t  +H) 
Z3-RV)BTO"  "AL3 


R(t  +H) 

Z4  -  C3F(tp.H)   -  C3  [1    -       ^tp)     3<AL4 

Goal   level    (initial   aspiration  levels) 


GL1   =  30.5  $/hr,  Gl_2  =  0.9885,           GL3  =  0.99,           GL4  =  $30,000 
Relevant  range  of  Z. 

r^)  '  [o,  80]  r(z2)  =  [0.9,  1] 

r(Z3)  =  [0.9,  1]  r(Z4)  =  [0,  200,000] 

Setup  procedure: 

We  begin  by  transforming  the  criterion  functions: 

Z,  Z9  -  0.9 

Y      x  —L  Y      =  - 

Tl       80  '  T2         1   -  0.9 

Z3  -  0.9  Z4 

Y„  =  -A =-=-  ,  Y,  -  - 


3  1   -  0.9   '  4       200,000 

The  initial   aspiration  levels  are  assumed  to  be  equal   to  the  goal   levels, 

AL.j   =  GL.j ,     i   =  1,2,3,4  and  so  the  values  of  A-   are: 

a      •  30-5  ..  n  iQioc  a      ■  0-9885  -  0.9      .  n  QQC- 

A-, g~ 0.38125  Ao ■,    _   r,  q =   0.885 


A     .  0.99   -  0.9  .     ■     30,000  _  Q  ]£ 

M3         1   -  0.9  u,y  rt4       200,000       U'i0 
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The  elements  of  d  are  formulated  from  the  structure  of  the  goals 


d]  =  Y^,        d2  =  A2/Y2,        d3  =  A3/Y3,        d4  =  Y4/A4 


First  Cycle: 


The  principal   problem  to  be  solved  on  the  first  cycle  is: 


min  S-j  =  d-|  +  d2  +  d3  +  d- 


s.t.         t    >  0 


Also  as  part  of  the  first  cycle,  we  construct  four  auxiliary  problems 
which  attempt  to  satisfy  each  of  the  goals  in  turn.  If  Z-,  <_  30.5  is 
entered  as  a  constraint,  d-,  is  deleted  from  the  surrogate  objective 
function  giving: 


min    s,  ,  =  d2  +  d3  +  d.  ' 


s.t.    Z1  <_  30.5 


t  >  0 
P 


Similarly,  the  second,  third  and  fourth  auxiliary  problems  are: 


min  s,   2  =  d-,  +  d3  +  d4 


s.t.        Z2  >_  0.9885 


*P  >0 
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min    s1  3  =  d-j  +  cL  +  d4 


s.t.    Z3  >  0.99 


t  >  0 
P 


min    sl  4  =  d]  +  d2  +  d3 


s.t.    Z4  <_  30000 


V  ° 


The  optimum  results  for  the  five  problems  are  tabulated  in  Table  4.2. 

The  information  from  these  five  problem  can  be  presented  to  the  DM 
in  the  graphical  form.  Figs  4.8  to  4.11  show  the  effects  of  imposing  goal 
attainment  at  the  current  aspiration  level  compared  to  solution  of  the 
principle  problem.  Let  us  compare  Z]  and  Z-, -,  in  Fig.  4.8.  The  shaded 
portion  of  each  axis  indicates  goal  attainment  relative  to  AL.  Assuming 
Zi  are  linearly  related,  then  dZg/dZj ,  dZ3/dZ] ,  dZ4/dZ]  are  constant. 
The  DM  uses  this  information  to  predict  the  approximate  levels  of  Z9,  I- 
and  Z4  for  a  given  aspiration  level  with  Z1  entered  as  a  constraint. 
Denote  a  desired  level  of  attainment  on  goal  1  by  Zy     Then  we  can  cal- 
culate the  effect  of  this  estimate  entered  as  a  constraint.  We  have: 

az2  =  (z1J2-z1<2)  .  (z1  -zM)/  (zK11  -  z,  ,) 


AV  (zi.i3-zi.3>  •  <zi  "zm)  /  (zi.n  -zM) 


53 


co 

O- 

o 

s — . 

s: 

o 

LU 

o 

CO 

o 

o 

CJ 

co 

!— 

u 

n 

5n 

o\ 

U 

ai 

■W 

o 

(/) 

s. 

« 

•1 — 

0.O 

4- 

CO 

co 

CJ 

en 

JC 

• 

4-> 

o 

<4- 

r> 

o 

m 

oo 

o 

+J 

P0 

r— 

3 

00 

ii 

aj 

eg 

o- 

CM 

ii 

* 

_i 

< 

CJ 

r— 

JO 

n 

Ml 


o 

oo 

CO 

O 

O 

CM 

<* 

en 

1X3 

o 

i — 

o 

-=3- 

CO 

o 

r 

ft 

r* 

#* 

ft 

r». 

o 

co 

co 

o 

i- 

on 

CO 

1 — 

CO 

a 

ft 

•V 

ft 

*N 

CNJ 

o 

co 

r^ 

o 

CO 

o 

«*■ 

CM 

co 

en 

co 

co 

en 

CO 

en 

en 

en 

en 

en 

• 

• 

• 

• 

• 

O 

O 

O 

o 

O 

'CO 

co 

on 

i — 

O 

on 

co 

co 

OO 

CO 

co 

co 

co 

co 

CO 

en 

en 

en 

en 

en 

■ 

• 

• 

• 

* 

o 

o 

O 

O 

O 

ft 

»> 

« 

ft 

ft 

CO 

en 

en 

CO 

a\ 

n-~ 

**■ 

M 

on 

OO 

• 

• 

• 

• 

• 

CM 

O 

CM 

. — 

=3- 

«■ 

CO 

CO 

«* 

oo 

II 

ll 

II 

ll 

II 

j 

•— . 

<N 

oo 

<*. 

r- ■ 

I-— 

^™ 

r^~ 

M~| 

M 

M 

M 

M 

r^ 

o 

O 

CNJ 

r>. 

^-^ 

*^N 

r^ 

- V 

on 

• — 

, — 

OO 

O 

■ 

■CO 

CO 

• 

> 

o 

OO 

o 

'CO 

CNJ 

C70 

r— 

r> 

• 

• 

ft 

ft 

CNJ 

1 — 

1 — 

CO 

en 

o 

1 — 

OJ 

o 

ft 

ft 

p- 

CNJ 

en 

OJ 

CO 

en 

O 

• 

on 

CO 

. 

• 

O 

CNJ 

oo 
O 

a 

r" 

ft 

. 

• 

n 

ft 

^O 

| 

^~ 

UO 

o 

r— 

r*. 

OO 

CO 

ft 

ft 

CNJ 

o 

o 

en 

o 

c 

o 

• 

oo 

• 

• 

• 

1 — 

en 

en 

"" 

■"" 

•"" 

ft 

• 

ft 

ft 

ft 

o 

O 

co 

on 

CO 

. — 

'CO 

r— 

r-. 

o 

ft 

on 

OO 

co 

<=r 

o 

O 

CO 

r— 

» 

II 

II 

II 

II 

II 

r— 1 

OJ 

en. 

■^", 

*~ ™ 

1 — ■ 

1— ' 

r-" 

r-» 

-a   I 

"O    1 

-o 

T3 

-o 

n 

co 

on 

on 

OO 

^— 

on 

OJ 

co 

r— 

OO 

O 

■en 

oo 

co 

U3 

cn 
0C 

en 

CO 

CO 

en 

-n 
O 

o 
o 

O 

CM 
OO 
KO 

i — 

CO 

CO 

CO 

CO 

CO 

'i 

ll 

ll 

OJ 

II 

CO 

II 

to 

c/> 

to 

CO 

00 

54 


200,000     1 


1       80 


O^E    0.9 


Fig.   4.8       First  cycle:     goal   1   satisfied  at  AL-|  =  30.5 
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200,000T      1   *      80T 


1*1 


0  £    0.91        0*         1° 
Z4         Z3         Zl  Y 


0.987- 


010.9 
Y2     Z2 


Fig.   4-9       First  cycle:     goal    2  satisfied  at  AL^  ■  0.9885 
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200,00^        %      80t         T  ] 


Fig.  4.10      First  cycle:     goal   3  satisfied  at  AI_3  =  0.99 
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1  T 200, 000 


Z14  Region  of  Acceptability 


minim  nit 


Fig.  4.11   First  cycle:  goal  4  satisfied  at  AL4  =  30,000 
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*V  (zi.i4-zi.4)  '  <zi  -zm}  I  (zi.n  -  zi-i} 

where  Z  ,  =  the  cth  element  of  Za  h: 
a. be  a.p 

Rather  than  ask  the  DM  to  solve  these  equations  for  trial  values  of  Z1 , 
we  can  solve  the  equations  graphically  by  Fig.  4.8.  Given  Z-j ,  we  ob- 
serve  that  lines  passing  through  itu  and  m4  find  the  Z^  and  Z2-  (The 
m2  and  m.  are  the  points  of  intersection  of  Z-j  and  Z-j  ^  in  Fig.  4.8). 
The  23  can  be  found  by  proportionality,  (i.e.  if  Z-j  is  the  midpoint  of 
(Z1  ^  -  Z,  ,),  then  Z3  will  be  the  midpoint  of  ( Z]  13-  Z1>3).  So  the 
DM  can  use  a  ruler  to  'try  out'  different  Z\  to  see  how  they  affect  the 
other  goals. 

Examining  these  results  shown  in  Table  4.2,  it  is  apparent  that  the 
change  in  availability  is  relatively  independent  of  attainment  or  non- 
attainment  of  the  other  goals.  It  seems  reasonable  to  choose  an  aspir- 
ation level  for  this  goal  and  enter  it  as  a  constraint.  From  Fig.  4-9, 
the  DM  assesses  the  impact  of  such  an  action  and  sets  a  new  aspiration 
level  of  0.987  for  goal  2.  DM  considers  that  such  a  decrease  is  not 
a  really  different  from  the  original  goal  of  0.9885. 

Second  Cycle: 

The  aspiration  levels  are  now  AL1  =  (30.5,  0.987,  0.99,  30,000) 
and  achievement  of  goal  2  has  been  entered  as  a  constraint.  The  principal 
problem  to  be  solved  on  this  cycle  is: 
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min  s2  2  =  d-,   +  cL  +  d. 


s.t.         Z2  >  0.987 


tp>  o 


and  the  three  auxiliary  problems  are: 


min  s2  ?!  =  d3  +  d4 


s.t.        Z2  >_  0.987 


Z1   <  30.5 


tp>0 


min  s2>23  =  d1  +  d4 


s.t.         Z2  >.  0.987 


Z3  >  0.99 


tp  >  o 
and 


min  s2>24  =  d1  +  d3 


s.t.         Z2  >_  0.987 


Z4  <  30000 


V° 
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The  optimum  results  are  presented  in  Table  4.3  and  their  effects  on 
goals  1,  3  and  4  are  illustrated  in  Figs.  4.12,  4.13  and  4.14. 

Inspection  of  Table  4.3  reveals  that  if  Z1  <  38.33  can  be  con- 
tented then  the  DM  will  reach  a  satisfactory  solution.  However  the  DM 
thinks  the  replacement  cost  per  unit  of  time  is  too  high.  From  Fig.  4.14, 
the  reduction  of  replacement  cost  can  be  achieved  by  increasing  the  ex- 
pected failure  cost  (goal  4).  We  assume  that  the  DM  has  modified  his 
aspiration  on  goal  4  as  Z,  <  50,000. 

Third  Cycle: 

The  aspiration  level   for  this  cycle  are  AL2  =  (  30.5,    0.987,  0.99, 
50,000)  and  both  goals  2  and  3  are  entered  as  constraints.     The  principal 
problem  is: 

min  s3>24  *  d1  +  d3 


s.t.         Z2  >  0.987 


Z„   <   50000 
4  - 


V° 


and  the  two  auxiliary  problems  are: 


mln  s3.243  ■  d3 


s.t.         Z2  >.  0.987 


Z4  ±  50000 


Z1  l  30.5 


V° 
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200,00Or      1  * 


Fig.   4.12     Second  cycle:     goal   1   satisfied  at  AL-.   =  30.5 
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200,000T       1«     80T         Tl 


.     ..  ©.. 


Oi:     0.9 

£     t: 


o  i 


It  1 


..   ©.. 


Fig.  4.13  Second  cycle:  goal  3  satisfied  at  AL3  =  0.99 
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1    T  200,000 


0.9       0.9 


Fig.  4.14     Second  cycle:     goal   4  satisfied  at  AL.f  =  30,000 
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and 


min  s3.243  =  d1 


s.t. 

Z2  >  0.987 

Z4  <  50000 

Z3  10.99 

t  >  0 
P 

The  results  are  tabulated  in  Table  4.4.     In  examining  Figs.  4.15  and  4.16, 
the  DM  learns  that  goal   1   and  goal   3  cannot  satisfied  simultaneously  at 
present  aspiration  level.     From  Fig  4.16,  the  DM  can  make  a  trade-off 
between  the  replacement  cost  per  unit  of  time  (goal   1)  and  the  mission 
reliability  (goal   3).     Assume  that  the  DM  has  modified  his  aspiration 
on  goal   3  so  that  he  will   be  contented  if  Z3  _>  0.985.     He  arrived  at 
this  decision  by  making  the  trade-off  between  goal   1   and  goal   3  on 
Fig.  4.16  and  can  roughly  predict  that  the  value  of  Z-|   for  Z3  =  0.985  is 
33.2  $/hr. 

Fourth  Cycle  and  Termination: 

The  principal   problem  to  be  solved  on  this  cycle  is 
min 

s.t. 


s4. 

,243 

=  dl 

h 

>  0. 

987 

Z4 

<_  500000 

h 

>  0. 

987 

s 

>  0 
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200,00Or        It      1  f 


Z 


1   T80 


1°-91 


Fig.  4.15  Third  cycle:  goal  1  satisfied  at  AL1  ■  30.5 
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200,000t      1   * 


1    *    1 


Fig.   4.16    Third  cycle:     goal   3  satisfied  at  AL3  =  0.99 
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and  there  is  no  auxiliary  problem.     The  results  are  tabulated  in  Table  4.5. 
From  this  Table,  the  DM  is  content  with  this  result  with  tp  =  912.9  hrs. 

Z1   =  C(tp)  =  32.52  $/hr 
Z2  =  A(t  )  =  0.9885 

Z3  =  R(tpH)   =  °'985 

Z4  =  C3F(tp,H)   =  $37,470 

4.5     Discussions 

In  section  4.1,  the  strictest  selection  method  was  applied  to  the 
sample  problem.     3y  selecting  the  shortest  preventive  replacement  age, 
the  minimal   requirement  values  set  for  some  criteria  were  all  satisfied. 
However,  if  there  were  the  cases  that  maximal  values  were  set  to  be  not 
exceeded  by  some  criteria,  then  the  strictest  selection  method  would 
mean  selecting  the  longest  preventive  replacement  age.     It  is  the  inten- 
tion of  the  strictest  selection  method  to  assure  that  all   the  minimal 
(or  maximal)  values  set  for  some  criteria  are  satisfied.     In  the  case, 
where  the  minimal   requirement  values  are  set  for  some  criteria  is  mixed 
with  the  maximal  values  are  set  for  other  criteria  in  a  problem,  the 
strictest  selection  will   be  infeasible.     The  D.M.   is  least  disturbed 
by  the  analyst  in  the  strictest  selection  method.     Once  the  D.M.  decides 
the  minimal    (or  maximal)   values   for  some  criteria,  the  solution  is 
essentially  determined.     Change  the  requirement  values  (minimal   or  maxi- 
mal)  for  some  criteria  is  the  only  way  the  D.M.   can  affect  the  solution 
of  the  problem. 
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For  the  lexicographic  method,  the  results  are  influenced  by  the 
DM's  ranking  of  the  importance  of  each  criterion.  The  solution  is  very 
sensitive  to  the  ranking  of  each  criterion.  Some  requirement  values 
(minimal  or  maximal)  for  lower  ranking  criteria  are  usually  not  satisfied. 

For  the  DM's  point  of  view,  the  Waltz's  lexicographic  method  is  very 
different  from  the  above  two  methods.  There  is  a  feedback  of  information 
to  the  D.M.  for  determining  the  allowable  relaxation  at  a  specific  value 
of  each  criterion.  Satisfying  solution  instead  of  optimal  solution  of 
each  criterion  is  pursued,  and  the  sensitivity  of  the  ranking  to  the 
solution  is  reduced. 

The  optimization  technique  used  in  SEMOPS  does  not  solve  problems; 
it  generates  information  so  that  the  DM  can  detect  the  inconsistent  con- 
straint set  and  select  an  acceptable  alternative.  The  key  feature  of 
the  SEMOPS  is  its  interactive  nature.  This  avoids  the  problem  of  speci- 
fying the  DM's  preference  structure  by  allowing  him  to  keep  within  himself 
the  transformations  that  he  makes  to  convert  numbers  into  value  judgments. 
The  DM  can  develop  a  ranking  of  goals  as  he  receives  information  concern- 
ing the  feasible  alternatives.  He  may  also  revise  his  preferences  during 
the  course  of  the  interactions.  For  the  numerical  example,  it  also  shows 
that  SEMOPS  prompts  the  modification  of  conflicting  aspirations  so  that 
an  acceptable  solution  can  be  determined.  SEMOPS  accomplishes  this  by 
revealing  to  the  DM  the  extent  to  which  his  aspirations  will  have  to  be 
modified  to  achieve  a  feasible  alternative. 
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These  four  methods  show  the  different  degree  of  participation  of  DM 
in  solving  the  multiple  criteria  decision  making  problem.  The  more  par- 
ticipation of  the  DM,  the  more  time  needed  in  analyzing.  Rather  to  find 
which  one  is  the  best  of  the  four  methods,  it  is  better  to  equip  with 
them  all  and  to  apply  them  in  following  manner: 

Solve  the  problem  by  the  strictest  selection  method  first.  If  the  DM 
satisfies  with  the  solution  obtained,  the  whole  problem  is  solved.  Other- 
wise the  analyst  should  ask  the  DM  to  supply  more  information,  such  as 
the  ranking  of  the  importance  of  each  criterion,  and  solve  it  by  the 
lexicographic  method.  If  the  DM  still  not  satisfies  with  the  second 
solution,  the  DM  should  be  asked  to  participate  more  in  the  decision 
making  process  for  determining  the  allowable  relaxation  at  a  specific 
value  of  each  criterion.  If  the  third  solution  is  still  considered  not 
satisfying,  then  the  interactive  method,  SEMOPS,  is  employed  to  generate 
information  to  the  DM.  By  which  he  can  detect  the  inconsistant  con- 
straint set  and  modify  the  conflicting  aspirations  so  that  an  acceptable 
solution  can  be  determined. 

It  becomes  obvious,  whether  the  more  complex  method  should  be  se- 
lected to  solve  the  problem  depend  on  the  DM's  satisfaction  with  the 
solutions  found  by  simpler  method  and  the  information  the  DM  can  supply. 
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CHAPTER  5.     CONCLUSIONS  AND  EXTENSIONS  FOR  FUTURE  STUDY 

This  study  is  to  determine  the  age  at  which  an  operating  critical 
item  should  be  replaced  with  a  new  item.     The  replacement  policy  is  to 
perform  a  replacement  once  the  critical   item  has  reached  a  specified  age 
(preventive  replacement)  or  a  failure  occurs   (corrective  replacement). 
A  variety  of  criteria  has  been  used  for  finding  the  optimal    replacement 
age.     However,   in  the  conventional   approaches,   usually  only  one  criterion 
was  considered.     It  is  obvious  that  the  replacement  age  will   depend  on 
what  criterion  is  used  to  evaluate  the  decision.     The  optimal   replacement 
age  for  one  criterion  is  surely  not  a  optimal   solution  for  other  criterion. 

This  study  employes  several  method  for  multiple  objective  decision 
making  to  determine  the  replacement  age  of  the  critical   item  of  a  main- 
tained system.     The  four  critera  considered  are:     (1)   the  replacement  cost 
per  unit  time,   (2)  the  availability  of  the  critical   item,   (3)  the  mission 
reliability,  and  (4)  the  expected  cost  of  failure  during  the  mission. 
The  solutions  are  found  by  using  the  following  four  methods  for  MO DM 
problems: 

(1)  strictest  selection  method 

(2)  lexicographic  method 

(3)  Waltz's  lexicographic  method 

(4)  sequential   multiple  objective  problem  solving  techniques 
(SEMOPS) 

Each  method  has  been  described  and  demonstrated  in  detail.     The  results 
obtained  by  each  method  are  directly  related  to  the  assumptions  and 
implications  of  each  method. 
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The  example  developed  in  the  present  study  is  an  one  decision  vari- 
able with  multiple  criteria  problem.  The  multiple  criteria  problem  with 
multiple  decision  variables  can  be  approached  by  the  same  way.  In 
addition,  the  models  for  the  replacement  of  a  single  item  can  be  extended 
into  group  replacement  models.  They  are  presented  below  for  the  future 
study. 

1 .  Multiple  Criteria  Problem  with  Multiple  Decision  Variables. 

For  the  problems  studied,  the  only  decision  variable  is  the  pre- 
ventive replacement  age,  t  .  The  cost  of  a  preventive  replacement,  C-j , 
the  cost  of  a  failure  replacement,  C-,  the  mean  replacement  time  for 
the  preventive  replacement,  t-.  and  the  mean  replacement  time  for  the 
failure  replacement,  t2>  are  all  assumed  to  be  independent  and  constant. 
However,  in  a  real  life,  the  replacement  cost  is  always  related  to  the 
mean  replacement  time.  If  shorter  mean  replacement  time  is  desired,  it 
requires  to  allocate  more  personnel  and  facilities  to  the  replacement 
work.  Therefore,  with  a  shorter  mean  replacement  time,  it  requires  a 
high  replacement  cost.  That  is,  the  costs  are  functions  of  replacement 
time  such  as 

C1  =  FUj)  (5.1) 

C2  =  F(t2)  (5.2) 

With  the  modification  of  (5.1)   and  (5.2),   the  sample  problem  in  this 
study  becomes  a  multiple  criteria  problem  with  three  decision  variables  : 
t  ,   t,   and  t«.     All   the  four  methods  demonstrated  in  this  study  is  appli- 
cable to  the  multiple  decision       variables  problem.     The  decision  process 
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is  unchanged.     Only  the  optimization  technique  should  use  optimization 
methods  for  multiple  variables. 

2.     Scheduled  Group  Replacement  Maintenance  Policy  Based  on  Multiple 
Objective  Decision  Making 

For  e^ery  piece  of  equipment  such  as  aircraft,  there  exists  critical 
item  whose  defect  could  result  in  hazardous  or  unsafe  conditions  for  in- 
dividuals using  the  equipment  or  could  prevent  the  successful   accomplish- 
ment of  the  equipment's  tactical   function.     As  a  consequence,   such 
critical   item  needs  to  be  replaced  before  the  failure. 

In  the  previous  study,   the  replacement  age  of  a  single  critical 
item  has  been  determined  subject  to  multiple  criteria.     Several  multiple 
objective  decision  making  methods  were  applied  to  determine  the  optimal   pre- 
ventive replacement  age.     The  solutions  and  implications  of  the  different 
multiple  objective  decision  making  methods  were  compared  and  explained 
in  detail . 

It  is  often  worthwhile  to  replace  similar  items   in  groups  rather  than 
singly,  since  the  cost  of  replacing  an  item  under  group  replacement  con- 
ditions may  be  lower,   i.e.   there  are  economies  of  scale.     The  determin- 
ation of  the  group  replacement  age  for  the  similar  critical   items  or  a 
group  of  non-similar  critical   items  should  be  also  based  on  multiple  cri- 
teria. 

The  proposed  problem  is   to  determine  the  optimal   group  preventive 
replacement  age  for  a  critical   item  of  a  equipment  based  on  the  following 
multiple  criteria: 


76 


(1)  the  replacement  cost  per  unit  of  time, 

(2)  the  availability  of  the  equipment, 

(3)  the  mission  reliability  of  the  equipment,  and 

(4)  the  cost  of  failure  during  mission. 

Assuming  there  are  several  similiar  items  which  are  subject  to 
failure.  Whenever  an  item  fails  it  is  replaced  by  a  new  item.  There 
is  also  the  possibility  that  the  group  replacement  can  be  performed  at 
fixed  interval  of  time.  The  item  replaced  will  behave  like  a  new  item 
with  zero  length  of  life.  Then  the  above  four  criteria  can  be  evaluated 
as  follows: 
(1 )  The  replacement  cost  per  unit  of  time 

If  f(t)  denotes  the  failure  density  function  of  an  item,  then  the 
reliability,  R(x),  and  the  unreliability,  F(x),  of  an  item  at  age  x  are 
respectively 

R(x)  =  /  f(t)dt 
x 

F(x)  =  1  -  R(x) 

Let  C  is  the  cost  of  replacing  one  item  under  conditions  of  group 
replacement,  Cf  is  the  cost  of  a  failure  replacement,  and  K  is  the  total 
number  of  items  in  the  group,  then  the  replacement  cost  per  unit  of  time 


C(tp)  is 


KC  +  KM  tn  C. 

C(t  )  =  — ^ P  f 

^  V     t„  +  Kt 
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where    t  =  the  group  replacement  age 

t  =  the  replacement  time  of  one  item  under  conditions 
9 

of  group  replacement 
M(t  )  =  the  expected  number  of  times  one  item  fails  in  interval 
(0,  tp). 

M(t)  can  be  determined  by  renewal  theory  [8]  as  follows. 

Let   N(t)        =  number  of  failures  in  interval  (0,t) 
t-.,t2,  ...,  t  =  intervals  between  failures 

S  =  time  up  to  the  rth  failure  =  r-|  +  r,,  +  ...  +  rp 

Fp(t)  =  cumulative  distribution  function  of  Sr- 
Now  the  probability  of  N(t)  =  r  is  the  probability  that  t  lies  between  the 
rth  and  the  (r+l)th  failure.  This  is  obtained  as  follows: 


P[N(t)  <  r]  =  1  -  Fp(t) 


P[N(t)  >  r]  =  Fr+1(t) 

since 

P[N(t)  <  r]  +  P[N(t)  =  r]  +  P[N(t)  >  r]  =  1 

P[N(t)  =  r]  =  Fr(t)  -  Fr+](t) 


M(t)  =  I     rP[N(t)  *  r]  =  J  r[F  (t)  -  F  ,(t)] 
r=0  r=0    r      r  ' 


M(t)  =  I     F  (t) 
r=l  r 

On  taking  Laplace  transforms  of  both  sides,  we  get 


7£ 


The  problem  is  then  to  determine  M(t)  from  M*(s).  This  can  be  done 
by  determining  f(t)  from  f*(s),  a  process  termed  inversion  Laplace  trans- 
forms. 

The  optimal  group  replacement  age,  t*,  that  minimizes  the  replace- 
ment cost  per  unit  of  time  is  obtained  by  solving 

dC(t.) 

—  P  =  0 
dt 
P 

(2)     Availability 

The  availability  of  the  equipment  when  all    the  critical   items  con- 
sidered are  operating,  A(t  ),  is  [13] 

t  -KM(t  )tf 

^y     tn  +  Kt 

p     g 

where  t  is  the  failure  replacement  time  of  one  item. 

The  optimal  group  replacement  age,  t*,  that  maximizes  the  availability 
is  obtained  by  solving 


dA(tn) 


*-■  0 


dt 
P 

If  there  is  requirement  for  a  standard  minimum  value,  A(t p)m-jn>   t0  De 
fulfilled  by  the  item,  then  it  is  required  that 

A(tJ    •     <  A(t*) 
v  p'min  -    v  p' 

and  the  highest  value  of  t*  satisfying  above  equation  should  be  selected. 
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(3)     Mission  reliability 

The  probability,  that  at  a  specified  time  T  the  unit  is  operating  and 
will  continue  to  operate  for  a  interval  of  time  H,  is  defined  as  "interval 
reliability  RI(H,X)"  which  is  given  by  [2] 

x 
RI(H,X)   =  R(X+H)  +  /     R(X-y+H)dMn(y) 

0 

where  dM,,(y)  is  the  probability  of  a  completed  repair  occurring  at  time  y. 

M..(t)  =  E[N..(t)],  the  expected  number  of  visits  to  state  j  in  (0,t)  if 

at  time  0  the  unit  enter  state  i.  The  operating  state  is  labeled  by  1  and 

the  failed  state  by  0. 

t 
Mn(t)  =  /  M01(t-z)dF(z) 


t 
0 


i01(t)  =  /  [1  +  Mn(t-z)]dG(z) 


where  F(t)  is  the  failure  time  distribution  and  G(t)  is  the  repair  time 
distribution.  If  F(t)  and  G(t)  are  known,  then  M]](t)  and  MQ1(t)  can  be 
determined.  The  mission  reliability,  R(H|X),  is  the  probability  thatall 
the  critical  items  are  operating  at  X  and  will  continue  to  operate  for  a 
interval  of  length  H. 

R(H|X)  =  R(H,X)K 

If  there  is  a  minimum  mission  reliability  requirement,   R(H|X)    .    , 

then  we  require  the  lowest  value,  i.e.,  R(H|t  )   still   higher  than 

R(H|X)mn-n 
v    '    'mm 
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R(H|X)m1„  <  R(H|t*) 

The  highest  value  of  t*  satisfying  the  above  inequality  should  be 
selected. 

(4)  Cost  of  failure  during  mission 

The  cost  of  failure  of  any  critical  item  during  the  mission,  C, 
includes  the  cost  of  damage  to  the  whole  equipment,  possible  loss  of 
lives,  etc.  It  has  a  higher  value  than  C  and  Cf.  The  expected  cost  of 
failure  during  the  mission  of  length  H,  E^,  is 

Ef  =  C[l  -  R(H|tp)] 

=  C[l  -  RI(H,tp)K]. 

The  condition  for  finding  the  optimal  group  replacement  age  t*   is  obtained 
from  the  requirement  that  the  expected  cost  of  failure  during  the  mission 
should  not  greater  than  the  group  replacement  cost  of  the  items  just  prior 
to  the  mission,  C  ,  plus  the  expected  cost  of  failure  of  the  newly  re- 
placed item  during  the  mission: 

C[l-R(H|tp)]  =  C[1-R(H|0)]  +  Cg 

Since  each  criterion  is  based  on  a  different  objective  function  and 
considers  different  factors,  this  problem  should  be  solved  by  methods  for 
multiple  criteria  decision  making  [7]. 
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Usually,  several  criteria  instead  of  one  should  be  considered  so  as 
to  determine  the  replacement  age  of  a  critical  item  of  a  maintained  system. 
In  this  study,  mathematical  models  have  been  developed  for  four  criteria: 
(1)  the  replacement  cost  per  unit  of  time,  (2)  the  availability  of  the 
critical  item,  (3)  the  mission  reliability,  and  (4)  the  expected  cost  of 
failure  during  the  mission.  The  solutions  are  obtained  using  four  methods 
for  multiple  objective  decision  making:  (1)  the  strictest  selection 
method,  (2)  the  lexicographic  method,  (3)  the  Waltz's  lexicographic  method, 
and  (4)  the  sequential  multiple  objective  problem  solving  technique 
(SEMOPS). 

Using  the  aircraft  engine  as  an  example,  the  optimal  replacement  age 
has  been  found  by  the  four  different  methods.  The  optimum  results  and  the 
implication  of  the  methods  are  discussed.  The  extensions  of  the  present 
study  to  the  multiple  criteria  problem  with  multiple  decision  variables 
and  to  the  group  replacement  problem  are  proposed. 


